High temperature Bose-Einstein condensation into an excited state at equilibrium 
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We describe Bose-Einstein condensation of strongly interacting particles into a quantum state 
which is an excited single-particle state, but becomes the many-body ground state as density in- 
creases because it minimizes the interaction energy compared to other states. Mean field calculations 
for a graphene potential just wide enough for two closely interacting layers of molecular hydrogen 
show condensation at temperatures up to 60 K. In the condensed state, molecules hop between 
layers, increasing the first peak in the pair-correlation function just past the hard core repulsion 
diameter. 
PACS numbers: 67.85. Hj, 67.85. Jk, 05.30.-d, 67.63.Cd 



A Bose-Einstein condensate (BEC) is a delicate state 
of ultra-cold matter in which a significant fraction of bo- 
son particles occupy the single lowest energy quantum 
state. While condensation was long suspected to play a 
role in superfluid helium [T] , the existance of a BEC was 
experimentally demonstrated [2jSj only recently. Since 
then, the primary application of the phenomenon has 
been in providing a macroscopic quantum object consist- 
ing of many particles acting in perfect unison, which can 
be probed in search of first principles for quantum me- 
chanics. A BEC only forms when the density of available 
quantum states (DOS) is reduced to the order of the ac- 
tual density, requiring temperatures T < 2.2 Kelvin (K). 

In this letter, we describe a new breed of BEC which 
forms from strongly interacting bosons that condense into 
a state with non-zero kinetic energy. The state mini- 
mizes the interaction energy, creating an energetic pref- 
erence for the condensate which is able to overcome the 
entropic cost of multiple occupation at higher tempera- 
tures than a traditional BEC. We show calculations using 
mean field theory for molecular hydrogen (H2) confined 
in a graphene potential with a precisely tuned width in 
which H2 forms two closely interacting layers. For an 
excited state in which a molecule hops between layers 
(FigflVc)), the first peak in the pair-correlation function 
[S] for two condensed H2 is increased past the hard core 
repulsion diameter [6] , reducing their interaction relative 
to any other pair of occupied states and producing a con- 
densate at temperatures up to 60 K. 

These relatively high temperatures will allow experi- 
mental studies under conditions which are considerably 
less delicate than those currently required. In addition, 
to our knowledge this is the first system for which bosons 
can be shown to condense into an excited state at equi- 
librium, without being driven by external conditions, as 
occurs with lasers [7; , condensates of quasi-particle exci- 
tations min], and vortices PMH]. solitons [13 E], and 
phonon-like collective oscillations [TS] evoked in typical 
BECs. As such, this condensate may provide a venue for 
new discoveries about quantum mechanics. 

Hydrogen adsorption potentials created by two paral- 
lel graphene layers^lG with AA stacking separated by 
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FIG. 1. Quantum states of H2 in a 8.1 A wide graphene 
slit-pore. The V=0 potential energy surfaces enclosing the 
adsorption volume are shown. Shaded ovals show regions of 
most concentrated probability density for three states: (a) a 
thermally averaged state at T =20 K (Color differences aid 
in distinguishing peaks.) (b) an aligned commensurate (AC) 
state (c) a staggered commensurate (SC) state. Remaining 
panels depict formation of SC state in Eq. ^ (d) A wavefunc- 
tion at the Dirac point. The magnitude at each hexagon cen- 
ter of the graphene is alteratively -1-1, 0, and -1. Alternation 
between the top and bottom occurs by combination of (e) the 
first two vibrational wavef unctions, (f) Adding/subtracting 
them concentrates the wavefunction at the top/bottom. 
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FIG. 2. Properties of H2 in a thermal distribution of states, or condensed in the aligned (AC) or staggered (SC) commensurate 
states. Pair correlation functions for two H2, at coverages of (a) c =0% and (b) c =67 % of a full monolayer of adsorbed H2 
for a w =8.1 A pore. The location of peaks in the pair-correlation function of bulk liquid H2O are also marked. Energy per 
particle as a function of pore width at (c) c =0% and (d) c =67 %. For the thermal distribution, the ground state energy is 
given. Energy per particle for a w =8.1 A pore as a function of (e) coverage and (f) condensate fraction in the staggered state. 



width w are highly corrugated in the potential energy 
along the plane, by ^45 K. As a result, the low-energy 
single-particle quantum states are heavily perturbed from 
their free equivalents. These states have peaks in their 
probability densities which form a triangular lattice with 
spacing a=2.46 A. In pores with w >7 A, the probabil- 
ity density also separates into two layers at the top and 
bottom of the pore, as depicted in Figure [ija) . 

In considering likely forms for the many-body quan- 
tum states, we rely on neutron diffraction experiments 
[TTJ HHj which show liquid II2 on graphene at T = 14- 
34 K forms an ordered phase commensurate with the 
underlying lattice, meaning the strong influence of the 
corrugation persists at high density. The lattice spacing 
a is smaller than the II2 hard sphere diameter, d=2.95 
A; below this separation the II2-II2 interaction becomes 
highly repulsive. Thus only every third lattice point is 
occupied, with an in-plane nearest neighbor distance of 
a' =4.26 A. The reciprocal lattice vectors of this trian- 
gular H2 lattice are equal to the wave-vector at the K 
and K' Dirac points of the underlying graphene lattice. 
In w ~7-10 A pores, the interaction between the top 
and bottom layers is also significant, and it is important 
to distinguish whether the two lattices are aligned (Fig. 
[ijb)) or staggered (Fig. ^c)). 

Matching quantum states may be easily formed from 
a few single particle states. We used the Mattera et al. 
[IB] C-H2 interaction to find the potential Vext{x) for 
a graphene pore, and solved the Schrodinger equation 
numerically |19) for the H2 center of mass using peri- 
odic boundary conditions for surface area A — 32 A^. 
The resultant quantum states are V'^ ^ n' "^here k gives 
the wavevector parallel to the plane in the first BriouUin 
Zone, ?Ti=l,2. . . is the band number, and n=0,l. . . is 
the quantum number of the confined vibrational motion 



perpendicular to the plane. Consider the following prob- 
ability densities, 
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where (3 — l/ki,T. These are the thermally averaged, 
aligned commensurate (AC), and staggered commensu- 
rate (SC) states depicted in Figure [T|a)-(c). Figure [TJd) 
shows the first two terms of Eq. |3] They describe a real 
sinusoidal wavefunction whose magnitude alternates be- 
tween -fl, 0, and -1 at adjacent lattice points, which is 
also in the ground vibrational state. The last term in Eq. 
[3] is the first excited vibrational state. Because the first 
two vibrational states (Fig. [lie)) can be combined to lo- 
calize the wavefunction at the top or bottom of the pore 
(Fig. [ijf)), the sum result is to create a wavefunction 
whose probability density hops between the top and bot- 
tom of the pore. The AC state is similarly formed, but 
without the combination of different vibrational states 
and thus without the hopping. 

As Figure [l] shows, adjacent peaks in the probability 
density have the largest separation for the SC state. In 
narrow pores where the distance between the two layers 
{I « w-b.S A) is small, the SC state is the only one whose 
adjacent peaks are further apart than d (Figl2ja)). As a 
result, the interaction of two II2 in this state is slightly 
smaller than for any other pair of states. 

When a high density of H2 is condensed into the SC 
state, the coherence between the spatial extent of the 
H2-II2 interaction and the natural periodicity of the SC 
state enhances the latter, futher reducing the interaction 



TABLE I. Parameters of the H2-H2 interaction functional. 
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energy for two H2 in this state. 

We calculated the many body quantum states using 
the mean field approximation, neglecting explicit cor- 
relation and exchange effects, and finding system ener- 
gies E{N,Na,Ns) and states as a function of number 
of total particles N, AC particles Na, and SC particles 
Ng, with the number of particles thermally distributed 
Nt — N — Na — Ng ■ These equations were solved itera- 
tively until self-consistent solutions were found. 
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Small changes in the numerical factors in front of each 
term in Eqs. [2][3] ensured the probability densities re- 
tained the same form. Numerically determined |19j so- 
lutions of Eq. [7] are repeated periodically for A' — 15A 
before use in Eqs. [4][5] Interactions are included through 
M, a finite-width H2-H2 interaction functional, 

(10) 
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similar to those used for liquid ^He [201 E] ■ The second 
term is increasingly repulsive with density and includes 
correlation effects in an approximate way. Parameters 
are given in Table llj e and a are typical Lennard-Jones 
parameters, [5] and the remainder reproduce properties of 
bulk liquid H2 [2ll[22]. 

We define coverage as c = 10.7A^/(yl x 2), where 10.7 
A^ is the surface area per H2 at monolayer coverage [T7] 
and the factor 2 accounts for both pore walls. Maximum 
coverages studied were the ones which hold exactly one 



H2 in each peak of the probability density of the SC and 
AC states, as long as all peaks are at least d apart (re- 
sulting in different cutoffs for each state as a function 
of w, e.g. Fig. l2]Je)). Above this density, it is certain 
that two H2 are often closer than d and we expect signif- 
icant correlation effects. The maximum c at any width 
was c = 67 %, equivalent to one commensurate lattice on 
each side of the pore wall. 

The pair correlation function for two particles in the 
condensate at this maximum coverage is shown in Figure 
^h). It is quite similar to that of bulk liquid H2.[5 (The 
major exception is the behavior at r < 2 A; this differ- 
ence originates from the mean field approximation.) As 
a result, the interaction of two H2 in this state is much 
smaller than for any other pair of states. 

Even though the single-particle SC state is always 
higher in energy than the single-particle ground state 
(Fig. [2](c)), the tendency of the SC state to reduce the in- 
teraction energy means that a condensate in this state be- 
comes the many-body ground state as coverage increases 
(Fig. [2](d)-(f)) for all pore widths studied. But it is only 
within a narrower width range, centered around w ~ 8 
A, that the energy difference between the SC state and 
a thermal distribution is large enough (Fig. W[d)) for a 
condensate to form at temperatures above the melting 
point of H2. It is in this range that the intermolecular 
interaction and the SC state are closely enough in phase 
to reduce the interaction energy significantly. 

Using the calculated energies E{N,Na,Ns) and the 
number of available quantum states. 



9t 



N,Na,Ns 



(T) = J2 exp[-P{El 

k.7n,n 



k.rn.n 



E, 



0,1,0 



)], 

(12) 

the canonical partition function was computed using 
Bose statistics, and interpolating between values of 
N,Na,Ns. The condensate fraction (/ =< N^ > /N), 
average system energy, and specific heat was computed 
from the partition function. The temperature (through 
j3 = 1/fcfcT) was fixed at T = 20 K in Eq. [l]and allowed 
to vary in Eq. 12 As a check, we performed calculations 
with T variable through the entire calculation for w = 8.1 
A, c = 67%, finding over T =20-80 K, only a <8% vari- 
ation in the relevant value E{Nt = iV) - E{Ns = N). 

Figure [3] shows the condensate fraction /, average en- 
ergy per particle, and specific heat as a function of T for 
several pore widths at a coverage of c = 67 %. There 
are similarities and differences between this system and 
a typical BEC. Because condensate formation relies on 
differences in the interaction energy, there is no sim- 
ple definition for the critical temperature Tc{w) based 
on the DOS. Nevertheless, the specific heat (Fig. IsFc)) 
has the same characteristic cusp as a BEC as / begins 
to grow, and falls toward the constant of a classical two- 
dimensional gas at T > Tc{w). 

If we define Tc{w) as the location of the cusp, we notice 
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FIG. 3. (a) Condensate fraction, (b) average energy per par- 
ticle, and (c) specific heat per particle as a function of tem- 
perature at 67% coverage for several different pore widths. 



Tc(7.3) = 60K is much greater than Tc{8.l) = 30 K, even 
though the energy preference for the SC state is larger in 
an 8.1 A pore (Fig. l2rd)). This is because the 7.3 A pore, 
having a smaller volume, has a smaller number of states 
gt ■ Thus the onset of condensation involves a balance be- 
tween the areal density of states, gt/A, which defines the 
entropic cost of multiple occupation, and the coherence 
of the SC state with the interparticle interaction, which 
constrains the energetic benefit of condensation. 

While we have not shown /(c), we note that the critical 
temperature varies less with c (i.e. density) than it may 
for a typical BEC. Because interactions are an essential 
part of this new condensate, c must be at least large 
enough for the interaction energy to become significant 
enough to make the condensate the many-body ground 
state, no matter what the temperature. 

We also find that the condensate fraction never fully 
reaches 1. This is because the energy per particle is 
roughly constant above a certain value of /, as shown 
in Fig. I2lf). At this critical value (which varies with w), 
the mean field of the condensate pulls all other states 
into forms which likewise minimize the interaction en- 
ergy and there is no longer a significant energetic benefit 
to increasing /. However, we have not included exchange 
in our calculations, a factor which tends to favor conden- 
sation into a single state [53] . 

We note that no condensate forms for w >9.3 A. Be- 
cause the properties of these larger pores approach those 
of independent graphene sheets, this result is consistent 
with studies of liquid H2 on graphene/graphite which 
have found no evidence of condensation or superfluidity 
[23] except below the melting temperature in the pres- 
ence of impurities [25j . 

Here we have described a unique method of creating 



a condensate of bosons. While we have presented cal- 
culations for a condensate of H2 which forms in narrow 
graphene slit-pores, it is likely that other species and/or 
different pore types may be substituted, as long as the 
same basic concept is employed: the pore must be de- 
signed so that a small subset of quantum states have 
a periodicity which mimics the high density bulk form, 
thereby minimizing the interactions and creating a strong 
energetic preference for those states. 

This method has several significant practical advan- 
tages when it comes to experimental study of Bose sys- 
tems. Because formation of a condensate in this sys- 
tem is not singularly dependent on the DOS and thus 
the temperature, the critical temperatures are orders of 
magnitude higher than those of other BECs formed of 
real particles. Not only would this reduce the sophisti- 
cation and cost of the cooling techniques required, but it 
also permits higher densities by eliminating the need of 
super-cooling. Equilibration during cooling and mainte- 
nance of the condensate during measurements are sim- 
plified both by the higher temperatures and densities al- 
lowed, as well as the thermal contact with the underlying 
structure which provides the trapping potential. How- 
ever, the optical techniques often used to study BECs 
will be obscured by the confining pore. Fortunately hy- 
drogen is well suited for scattering studies due to its high 
neutron cross section, and deep-inelastic scattering has 
been developed as a technique to study BECs in helium, 
both in bulk form [35] as well as when adsorbed [?7] . 

While lasers with A wavelengths have been demon- 
strated [7,, it is likely that formation of the described 
pore with graphene will be more easily achievable than 
with an optical lattice. Previous work has demonstrated 
that graphene slit-pores with controllable widths may 
be formed from graphene oxide [55] sheets connected 
by benzene-l,4-diboronic acid linkers, with linker den- 
sity controlling width. Fortunately, the w ~ 8 A range 
requires a small density of linkers, and if oxides can be 
adequately cleaned from the surfaces, this system may 
represent a likely candidate in which the Bose conden- 
sate described here may be observed. 
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